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Glasses 1,2 are amorphous solids, in the sense that they display elastic be- 
haviour. In crystals, elasticity is associated with phonons, quantized sound-wave 
excitations. Phonon-like excitations exist also in glasses at very high frequencies 
(THz), and they remarkably persist into the supercooled liquid 3 . A universal 
feature of these amorphous systems is the Boson peak: the vibrational density 
of states g{oS) has an excess over the Debye (squared frequency) law, seen as a 
peak in g(uj)/uj 2 . We claim that this peak is the signature of a phase transition in 
the space of the stationary points of the energy, from a minima-dominated phase 
(with phonons) at low energy to a saddle-point dominated phase^»£ (without 
phonons). Here, by studying the spectra of inherent structures 7 (local minima 
of the potential energy), we show that this is the case in a realistic glass model: 
the Boson peak moves to lower frequencies on approaching the phonon-saddle 
transition and its height diverges at the critical point. The numerical results 
agree with Euclidean Random Matrix Theory 8 predictions on the existence of a 
sharp phase transition 9 between an amorphous elastic phase and a phonon-free 
one. 

Vibrational excitations of glasses in the THz region, and the related vibrational density 
of states (VDOS), play a crucial role in their thermodynamic properties^. Recent results 
suggest that the VDOS is determined by the properties of the potential energy landscape 7 
(PEL), a crucial tool to understand the slow dynamics of supercooled liquids 2 . Indeed, the 
VDOS and the dynamic structure factor can be qualitatively reproduced from the (harmonic) 
vibrational spectrum obtained from the diagonalization of the Hessian matrix of the potential 
energy evaluated at the minima (or inherent structures [IS]) of a Lennard- Jones system^. 
The same numerical procedure led to quantitative agreement with inelastic X-ray scattering 
experiments in amorphous silica (O. Pilla et al., preprint condmat/0209519). Here we will 
show that a Boson peak must be present in the VDOS as a consequence of the PEL geometry, 
and make quantitative predictions which can be checked fastly quenching supercooled liquids 
below their glass temperature. 

The high-frequency (0.1-10 THz) excitations have been experimentally shown to have 
linear dispersion relation o 3 i 12 i 13 i 14 i 15 i 16 i 17 in the mesoscopic momentum region (~ 1-10 nm _1 ). 
Although clearly not plane waves, they are highly reminiscent of phonons because they 
propagate with the speed of sound and because the VDOS, g(u), is Debye-like (g(uj) oc u 2 ) at 
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low enough frequency. These excitations have in fact been dubbed "high-frequency sound" 
and a coherent theoretical picture of their properties was obtaine d 9 ' 18 ' 19 using Euclidean 
Random Matrix Theory 8 (ERMT). 

However, there is more to the low frequency VDOS of glasses than the Debye law. 
A characteristic feature is that the VDOS departs from the Debye form, displaying 
an excess of states, which has been named "Boson peak" (BP). As observed in many 
material o 3 ! 12 ! 13 ! 14 ! 15 ! 16 ' 17 , the BP is in a region of frequencies where the dispersion relation for 
phonons is still linear (uo ~ a few THz). Several non-equivalent ways of experimentally defin- 
ing the BP have been introduced. One approach defines it as a peak in Raman-scattering 
data, another as a peak in the difference between the VDOS of the glass and the corre- 
sponding crystal. Other authors extract g(u>) from their neutron scattering data and look 
for a peak in g(u)/uj 2 . This is the definition we adopt, since we believe it unveils universal 
properties of glasses. 

The origin of the Boson peak can be understood if we consider the ensemble of general- 
ized inherent structures (GIS): for each equilibrium configuration the associated GIS is the 
nearest stationary point of the Hamiltonian. If we start from an equilibrium configuration 
at low temperature, the GIS is a local minimum, and it coincides with the more frequently 
used IS (i.e. the nearest minimum of the Hamiltonian). On the contrary, if we start from 
high temperature, the GISs are saddle points. In the GIS ensemble there is a sharp phase 
transition separating these two regimes. It takes place in glass-forming liquids*^ at the 
Mode Coupling temperatur e 21 ' 22 (Tmc), above which liquid diffusion is no longer ruled by 
rare "activated" jumps between ISs but by the motion along the unstable directions of sad- 
dles. Phonons are present in the spectrum of the VDOS in the low temperature phase (IS 
dominated) but are absent in the saddle phase (see ref. |23jfor a different approach). 

The key point, following both from analytic computations in soluble models 24 and 
simulations^ 1 *^, is that the minima obtained starting from configurations below Tmc and 
the saddles obtained starting above Tmc join smoothly at Tmc- Thus we can study GIS as 
a single ensemble, parametrized by the initial temperature or the final energy. We expect 
that the GIS ensemble belong to a large universality class. 

This transition from the phonon phase to the saddle phase is a quite general phenomenon 
that it is not restricted to glasses. It has been studied in the framework of the Euclidean 
Random Matrix Theory (ERMT)* 8 -, showing 9 that close to this transition, a BP is present in 
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the phonon phase, whose position shifts to lower frequencies on approaching the transition 
point, while its height grows without bound. The BP actually signals a crossover at frequency 
uj bp between a (phonon dominated) u 2 scaling of g(uS) to an u 1 scaling (7 < 2) that is 
present at the phase transition point. More precisely, at frequencies small respect to the 
Debye frequency, the VDOS should satisfy the scaling law 

g(u,A)=u<h(uA-P), (1) 

where A is the distance from the critical point and depends on the actual control parameter 
(pressure, temperature, etc.) and h(x) is such that h(x) ~ x 2-7 for x <C 1 and h(x) ~ 
constant for x ^> 1. This scaling law implies that 

-bp oc A p , oc A-?, (2) 

where (3 = p(2 — 7). Under the resummation of a given class of diagrams, ERMT predicts 9 
that p = 1, 7 = 3/2 and (3 = 1/2. Of course the actual (universal) values of these critical 
exponents in three dimension may differ slightly from the values found in this approximation. 

We now analyze numerically the existence of such a transition by studying the soft sphere 
binary mixture, a model of a fragile glasst^. We simulate this system with the Swap Monte 
Carlo algorithm 26 , and compute the VDOS of the ISs obtained starting from equilibrium 
configurations at temperatures below T MC . Except were stated, we used 2048 particles, and 
runs were followed until the energy of the IS, eis, reached a stationary value. The predictions 
above agree with the numerical data, taking A = e c — eis (eis being the energy of the ISs 
and e c the critical value). 

Given the numerical results of Kob et al?2£ one expects that g(u) depend only on eis, 
and thus that eis be the relevant control parameter. To check this, we compare the VDOS 
of two systems at different temperatures (one in equilibrium, the other not) but with the 
same e^. The equilibrated system has 2048 particles, while the other has 20000 and is in 
the regime where eis has not yet stabilized after a quench (Fig. lb). As Fig. la shows for 
a representative case, the VDOS obtained at the same eis f° r the two systems coincides, 
confirming that eis is the control parameter. We can thus proceed to investigate the scaling 
laws (j2J). 

The VDOS at three representative temperatures is shown in Fig. 2. For temperatures up 
to T = 0.69 T MC , a Debye behaviour (g(uj) oc u 2 ) is found for the lowest frequencies that can 



be resolved in a system of this size. In the g{uj)/uj 2 plot (Fig. 2a), a peak is clearly identified, 
which is seen to grow in height and shift to lower frequency on rising the temperature. The 
peak shifts to frequencies below our resolution for T > 0.69 T MC , but we can identify a Debye 
region up to T = 0.83 Tmo Above this temperature the Debye region is not well defined, and 
it disappears completely for T > 0.9 T MC . At these temperatures, g(u)/u 2 seems to diverge 
(at least within the frequency range that we can resolve). Experimentally, this behaviour 
may be difficult to observe due to the presence of the elastic peak, but B2O3 data2& seem 
to support this prediction. For frequencies immediately above cj B p, the VDOS scales as 
cu 7 , with 7 m 1.5 (Fig. 2b). Note that these changes of the VDOS with temperature go 
practically unnoticed if one plots g(oo) directly (Fig. 2c). In particular, the BP is completely 
unrelated to the first maximum of g(w), whose position is insensitive to temperature changes. 

Since we know that eis decreases slowly with time after a quench^! (see also Fig. lb), 
the above results enable us to make predictions about the ageing of the VDOS, and in 
particular of the BP. With increasing time the system moves farther from the critical point, 
and thus the BP should decrease in height and shift to higher frequencies. Moreover, at a 
given frequency (below the asymptotic cupp) g{w) should decrease, since it will be of order 
uj 1 at short times and of order uj 2 at very long times. Similarly, one should expect a cooling 
rate dependence of the shape of the BP: the slower the cooling, the lower the asymptotic 

and thus a larger uj-qv and less pronounced BP. An effect of this kind has been observed 
in AS3S3 28 , although here it is difficult to disentangle the physical effect from the chemical 
changes of the sample due to dissociation. To our knowledge, a clear-cut experimental 
observation of these effects has not been reported. 

Using all the spectra for which the peak position can be clearly identified, we find that the 
relationship between uj-qp and the energy of the IS is linear (Fig. 3a). The energy at which 
ubp becomes zero, e c , is found from a linear fit as e c = 1.74 e (e is the energy scale). As 
for the height of the peak, the points up to T = 0.83 Tmc are compatible with a power-law 
divergence of the form (J2J) (Fig. 3b). Fixing e c at the value 1.74 e arising from the linear fit of 
ubp vs. eis, a power-law fit yields an exponent (3 = 0.40(15), while if one fixes the exponent 
(3 = 1/2, then the critical value is obtained as e c = 1.752(2) e. Thus the numerical data are 
compatible with the theoretically predicted scaling, although we have not been able to work 
close to the critical point, and thus cannot get a great accuracy on the critical exponents or 
the critical point. 
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The present discussion applies to experiments as long as one is in the regime where the 
inverse frequency is much larger than the structural relaxation time, when the harmonic 
approximation makes sensed. Thus, while the BP observed in liquids can be understood 
by the present considerations, one should not expect to actually arrive close to the critical 
point in equilibrium measurements. Rather, one should generate a glass by hyperquenching 
(see e.g. talk by C. A. Angell in http://www.df.unipi.it/workshop/oral/OL6.pdf): for this 
we predict that the BP measured at low T should obey the scaling (j2J), the control parameter 
being the Active temperature, which in our instantaneous quench is just the temperature at 
which the hyperquench starts, but in experiments involves also the cooling rate. In this way 
one should also alleviate the problem of the swamping by the quasielastic peak. 

The significance of the value e c = 1.74(1) e is twofold. First, it is the value of eis that 
corresponds, in equilibrium, to T = Tmc, the Mode-Coupling critical temperature (which we 



have independently determined as in ref. 



22, using a standard Monte Carlo without swap). 



Second, it coincides within errors with the threshold energy below which the ratio of the 
number of saddle points to the number of minima vanishes exponentially with the number of 
particles^ (our numerical accuracy on e c is similar to that of|6J). Thus we find a link between 
the onset of slow dynamics, the appearance of phonon-like excitations and a geometrical 
phase transition of the PEL in fragile glasses. At variance with previous PEL studies, the 
clear link between the PEL geometrical transition and high-frequency dynamics allows it to 
be studied experimentally by following the evolution of the BP. 

In summary, we have studied numerically the VDOS of a fragile glass- former liquid under- 
going its Mode-Coupling transition. We have found that in the ensemble of the stationary 
points near the equilibrium configurations (GIS) there is a transition from a mechanically 
unstable phase where saddles are present to a stable phase where phonon-like excitations 
appear. The scaling laws predicted by ERMT are compatible with the numerical results we 
have obtained for the soft sphere model. This approach is equivalent to other microscopic 
descriptions of the glass transition based on the geometry of the potential energy landscape 
(PEL ) 4 '^'^ . But this point of view emphasizes quantities accessible to experiment, like the 
VDOS, and proposes an interpretation of an universal feature of glasses, the Boson peak. 
The BP appears in the phonon phase, and it is a signature of a cross-over from a phonon- 
dominated spectrum with a Debye ui 2 scaling to an w 7 , 7 » 1.5 spectrum, resulting from 
the hybridization of acoustic modes with high-energy modes that soften upon approaching 
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the saddle-phonon transition^. The ERMT predictions (eqs. ^ and |2J) could be checked ex- 
perimentally in hyper quenching experiments, where the Boson peak should strongly depend 
on the Active temperature. We expect that the saddle-phonon transition point of view will 
be able to bridge the realms of experiment and numerical studies of the PEL, allowing to 
test many geometrical ideas and to use the PEL insights in the detailed analysis of the 
experimental glass transition. 
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FIG. 1: The VDOS g(u>) at low frequencies depends only on the energy of the inherent structure, eis- 
a. g(uj)/uj 2 for an equilibrated 2048 particle system at T = 0.9 Tmc compared to that obtained 
from ten inherent structures of the same eis from a 20000-particle system out of equilibrium. 
The larger system was quenched from infinite temperature to T = 0.54 Tmc and the inherent 
structures obtained by minimizing instantaneous configurations between 1.9 • 10 4 and 2 • 10 4 Swap 
Monte Carlo steps. Similar results are obtained at all the temperatures for which comparable 
eis a r e obtained in the quench. IS were found by minimizing instantaneous configurations with 
a conjugate gradient algorithm. The Hessian was diagonalized with standard library routines in 
the case of the small system, and with the method of moments for the big one. Each VDOS was 
obtained from at least 200 IS as the histogram of the square root of the eigenvalues, b Monte 
Carlo history of eis i n logarithmic time-scale, for the system of 20000 particles, suddenly quenched 
from infinite temperature to T = 0.54 Tmc- The soft-sphere binary mixture is made of 50% of 
particles of type A and 50% of type B, both with the same mass. The interaction potential is 
Vaj3( r ) = + &f3)/ r ] 12 + («, P £ {A -£>})) wrtn cr's fixed by the conditions <jb = 

and (2cja) 3 + 2((ta + cb) 3 + (2<tb) 3 = 4<7o. A smooth cutoff is imposed at r c = \/3cro : f° r 
r c < r < a, we set V a p = B a p(a — r) 3 , and V a p = for r > a. In the results, lenghts are given 
in units of do, energies in units of e and frequencies in units of ujq = (mcjg/e) -1 / 2 . Using Argon 
parameters (gq = 3.4A, e = 120 K/cb, m = 39.96 u.a.m.) the frequency unit ujq is 0.46 Thz and 
the Mode Coupling temperature Tmc is 26AK. All runs are at constant volume, with particle 
density p = cr^ 3 . For temperatures down to T = 0.9 Tmc thermalization was achieved, as checked 
by ensuring that the equilibrium fluctuation-dissipation ratio is reached. For T < 0.9 TmC) the run 
was followed until eis reached stationarity. 
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FIG. 2: Vibrational density of states g(uS) of the soft-sphere binary mixture at three representa- 
tive temperatures (the full set of temperatures was T/Tmc = 0.9,0.83,0.78,0.69,0.61,0.54,0.49, 
frequencies are given in units of uj ). a. The g(uj)/uj 2 plot shows a Boson peak at low temperature 
(shown T = 0.54 Tmc)- When temperature is risen, the peak shifts to lower frequencies and grows 
without bound, b. The divergence seen in a is due to an w 7 (with 7 < 2) scaling of g(u>). At tem- 
peratures where the Boson peak is still seen, a crossover between J 2 and w 7 scaling can be noted 
comparing a and b. While 7 cannot be extracted very accurately from the data, it is compatible 
with 7 = 2 — (3/p w 1.6, with [3 and p taken from the fits of Fig. 3. c. The high frequency features 
of g(u>) do not change drastically with temperature. 
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FIG. 3: Scaling of the position, ujbp, and height of the Boson peak near the saddle-phonon transition 
(energies and frequencies in units of e and ujq respectively), a. wbp is linear in the control 
parameter, in this case the energy of the inherent structures eis- Extrapolating with a linear fit, 
ujbp goes to at eis = e c = 1.74(1) e. b. The height of the Boson peak diverges as a power law (the 
leftmost point has been left out of the fit). Height and position of the BP were obtained by fitting 
a parabola to the peak of g{uj)/uj 2 at T/T MC = 0.49, 0.54, 0.61, 0.69, 0.78, 0.83. The corresponding 
eis/e are 1.668, 1.671, 1.674, 1.682, 1.689, 1.707. 
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